的 变化 对 吸力 
马赫 数 为 0.4。 而 
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Abstract: In order to predict and analyze the phenomenon of separation, reattachment as well as separation 


bubble at different incidence angles. Self-developed large eddy simulation (LES) program of compressible 


N-S equations was used to simulate the flow field and the evolutionary process of the boundary layer in 


T106A low-pressure turbine cascade. The outlet Mach number is 0.4 and the Reynolds number is 1.1x10° 


based on the exit isentropic velocity and the axial chord length. The computational results show the 


time-averaged pressure coefficient on the blade surface and the points of separation and reattachment at 


+7.8° incidence angle compare well with experiment data. The boundary layer separates from the suction 


side and the hairpin vortexes are formed in succession. The rapid growth of the hairpin vortex eventually 


leads to transition in boundary layer. As the incidence angle changes from positive to negative, the separation 


point on the suction side moves downstream and the length of the separation bubble is gradually shorter. 
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国内 外 很 多 学 者 通过 大 量 的 实验 和 数值 模拟 广 
0 前 言 法 ， 分 析 了 攻 角 对 涡轮 叶 栅 流 场 结构 和 性 能 等 方面 
_ 影响 。Jouini 中 和 Benneri 等 人 都 通过 实验 发 现 ， 
头 EF [al 在 大 培 具 、 上 £ 包 

4 攻 角 在 很 大 程度 上 影响 着 叶 叶 型 损失 及 二 次 流 扣 
备 良好 的 工作 能 力 。 低 压 涡 轮作 为 航空 发 动机 的 主 。 失 的 大 小 。Dossena 等 ”的 研究 结果 表明 ， 气 流下 角 
要 部 件 之 一 常 工作 王 低 雷诺 数 人 的 变化 会 引起 二 次 流 大 小 、 位 置 和 形状 的 改变 。 张 
En | 获 等 四 对 雷诺 数 为 5x10 时 的 涡轮 叶 棚 进行 模拟 ， 描 
交角 是 影响 油轮 时 楼 变 工 况 气动 性 能 的 主要 因素 之 。 从 守 攻 m2 二 全 入 

一 ， 详 细 研 究 不 同 攻 角 下 低压 涡轮 叶 栅 流 场 的 非 定 。 全 从 ee 
值 模拟 对 二 维 涡轮 叶 栅 流 场 进行 研究 ， 结 果 显 示 攻 


常 特性 具有 重要 的 意义 。 
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角 对 近 壁 面 流动 的 影响 主要 体现 在 前 缘 附 近 。 陈 绍 
文中 ， 孙 大 伟 站 等 也 对 涡轮 流 场 、 损 失 及 性 能 的 变 攻 
角 特 性 进行 了 探讨 。 这 些 工 作 主要 集中 在 分 析 攻 角 
变化 对 通道 涡 、 二 次 流 损失 或 总 压 损 失 所 产生 的 影 
响 ， 而 对 流 场 结构 非 定常 演化 过 程 的 研究 较 少 。 

攻 角 的 改变 会 引起 分 离 泡 形态 和 尺寸 的 变化 ， 
进而 影响 涡轮 性 能 。 大 涡 模 拟 方法 凭借 其 独特 的 优 
势 ， 能 够 较 好 的 描述 流 场 中 流动 细节 的 非 定常 特性 。 
本 文采 用 自 编程 三 维 可 压缩 流动 的 大 涡 模 拟 方法 ， 
对 不 同 攻 角 下 涡轮 叶 栅 非 定常 流 场 进行 计算 ， 研 究 
其 对 边界 层 分 离 再 附 的 影响 ， 分 析 边 界 层 演化 过 程 ， 
为 涡轮 的 优化 设计 提供 更 多 理论 基础 。 


1 控制 方程 与 数值 方法 


本 文 所 用 的 控制 方程 和 数值 求解 方法 与 文献 [10] 
基本 一 致 ， 对 流通 量 通过 四 阶 精 度 的 中 心 格式 计算 ， 
亚 格 子 模型 采用 动力 涡 粘 模型 ， 时 间 推 进 为 三 阶 三 
步 的 Runge-Kutta 法 。 计算 程序 为 多 块 并 行 大 涡 模拟 
程序 (MPLES)， 其 更 多 细节 可 参见 文献 [10]。 模 拟 
+ 机 中 的 流动 属于 内 流 问 题 ， 边 界 条 件 通常 给 定 进 
口 总 温 、 总 压 、 气 流 角 以 及 出 口 背 压 ， 由 文献 [11] 给 
出 的 方法 求 得 。 


2 计算 模型 


低压 涡轮 叶 栅 T106A 的 叶 型 及 几何 参数 分 别 如 
图 1 和 表 1 所 示 。 图 2 为 计算 域 示意 图 和 边界 条 件 
的 设置 方式 ， 计 算 域 入 口 距 叶 栅 前 缘 0.5 Cj (Ci 为 
轴 向 弦 长 ) 出 口 距 叶 栅 尾 缘 1.6 C,,, 展 向 高 度 为 0.2 
Cuw。 无 量 纲 化 以 后 的 流向 、 周 向 及 展 向 尺寸 分 别 为 
3.1、0.93 和 0.2。 三 个 方向 的 网 格 为 320x180x60， 
壁面 第 一 层 网 格 y<1， 边 界 层 内 布置 约 35 个 节点 。 
基于 轴 向 弦 长 和 出 口 速度 的 雷诺 数 为 1.1x10”， 出 口 
马赫 数 为 0.4， 来 流 攻 角 的 变化 范围 是 -10° ~ +7.8°。 


七 


30.7 


axial chord 


图 1 T106A 叶 棚 叶 型 
Fig.1 Profile of T106A cascade 
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表 1T106A 叶 栅 几何 参数 


Table 1 Geometrical parameter of T106A cascade 


参数 单位 数值 
纺 长 [mm] 198 
轴 向 弦 长 [mm] 170 
安装 角 [9] 30.7 
节 距 [mml] 158 
叶 高 [mm] 375 


Inlet 


[人 SS 
。 
Suction Slide 


Periodic Boundary 


Periodic Boundary 


Outlet 


图 2 计算 域 示 意 


Fig.2 Computational domain 


3 计算 结果 分 析 

当 流 场 充 分 发 展 ， 即 流动 参数 具有 准 周 期 性 的 
变化 规律 时 ， 开 始 采集 瞬 态 流 场 ， 每 个 工 况 的 总 采 
样 瞬 态 流 场 均 多 于 3000 个 ， 总 采样 时 间 均 超过 26 
个 周期 。 对 叶片 沿 周 向 所 受 合力 随时 间 的 变化 进行 
快速 传 里 叶 变 换 ， 图 3 给 出 了 +7.8? 攻 角 下 叶片 沿 周 
向 受 力 的 频谱 图 。 图 中 可 以 看 出 , 主 频率 为 广 4.298， 
因此 无 量 纲 的 准 周期 为 7=0.2327。 其 他 来 流 攻 角 工 
况 的 主 频 率 和 准 周期 值 与 +7.8° 攻 角 结 果 相差 不 大 ， 
在 表 2 中 给 出 。 


表 2 各 攻 角 工 况 的 主 频率 和 准 周期 


Table 2 Main frequency and period of different incidences 


攻 角 主 频率 (7) 准 周期 (7) 
-10° 4.549 0.2198 
-5° 4.733 0.2113 
0° 4.342 0.2303 
+5° 4.442 0.2251 
+7.8° 4.298 0.2327 


图 4 为 +7.8° 攻 角 下 叶 栅 表面 时 均 压 力 系数 与 实 
验 测量 数据 的 对 比 ， 压 力 系 数 的 定义 为 : 


cys=(p-p,)/(p; —p;) (1) 


式 中 pi “和 ps 分 别 为 入 口 总 压 和 出 口 静 压 。 可 以 看 出 ， 
本 文大 涡 模拟 结果 与 实验 数据 吻合 较 好 ， 只 是 对 吸 
力 面 后 部 压力 平台 的 预测 略 有 区 别 ， 可 能 是 由 于 数 
值 模拟 与 实验 的 来 流 满 流 强 度 有 所 差异 造成 的 。 吸 
力 面 后 部 的 压力 平台 说 明 ， 壁 面 附近 流体 向 下 游 输 
运 过 程 中 ， 受 逆 压 梯度 和 粘性 力 等 因素 的 共同 作用 ， 
流速 逐渐 降低 ， 发 生 了 边界 层 分 离 。 

从 旋涡 角度 分 析 边 界 层 的 分 离 、 转 换 及 再 附 等 
现象 ， 有 助 于 了 解 边界 层 的 发 展 规律 ， 图 5 给 出 了 
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4.298 
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3 +7.8° 攻 角 下 叶片 受 力 的 频谱 图 


Fig.3 FFT of the force on the cascade 


Experiment 


图 4 叶 栅 表面 时 均 压力 系数 分 布 


Fig.4 Time-averaged pressure coefficient distribution 
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+7.8? 攻 角 时 ， 某 两 个 时 刻 吸力 面 后 部 大 尺度 相干 结 
构 的 形态 ， 采 用 2 法 则 进行 旋涡 识别 ， 取 CO=120 的 
等 值 面 。 图 中 可 以 看 出 ， 上 游 发 展 来 的 二 维 剪 切 层 ， 
由 于 展 向 速度 扰动 的 影响 迅速 表现 出 三 维特 性 ， 展 
向 涡 发 生 扭 曲 ，4 时 刻 形成 了 三 角形 的 和 A 涡 结构 。 入 
涡 向 下 游 输 运 过 程 中 不 断 生 长 ， 很 快 便 呈 现 出 类 似 
于 发 卡 涡 的 形态 (如 图 Sb) 所 示 )。 发 卡 涡 的 拉 伸 作 
用 逐渐 导致 了 层 流 边 界 层 的 裔 溃 ， 大 尺度 的 旋涡 随 
即 开 始 破碎 ， 形 成 大 量 较 小 尺度 的 旋涡 ， 小 尺度 旋 
ss 边界 层 内 流体 发 生 转 
再 附 到 壁面 上 ， 流 动 完 全 变 为 清流 状态 。 


ee <€ 
& SS 


荡 , 


(DA (b) 亡 


图 5 吸力 面 后 部 的 大 尺度 逆序 结构 


Fig.5 Coherent structure at the rear part of suction side 


图 6 为 +7.8°? 攻 角 工 况 下 ， 吸 力 面 瞬时 切 应 力 的 
时 衬 分 布 图 ， 图 中 黑色 实 线 为 切 应 力 等 于 零 的 位 置 。 
可 以 清晰 地 看 到 ,在 ee 
分 离 ， 但 在 某 些 时 刻 ， 分 离 的 边界 层 会 立即 附着 到 
壁面 上 ， 随 后 再 次 分 离 ， 从 而 形成 了 多 个 分 离 泡 同 
时 存在 的 现象 。 

图 7 为 不 同 攻 角 下 吸力 面 时 均 摩 探 力 系数 的 分 
布 。 摩 探 系 数 的 定义 为 : 


C=7,/(p;-p;) (2) 


式 中 ps 为 出 口 总 压 ，t 为 沿 轴 向 的 壁面 切 应 力 。 首 
先 观察 +7.8° 攻 角 工 况 ， 图 中 显示 ， 从 前 缘起 摩擦 力 
系数 先 减 小 到 一 个 极 小 值 ， 随 后 增 大 ， 这 与 图 4 中 
吸力 面 压 力 系数 先 增 大 后 减 小 的 变化 相对 应 ， 叶 栅 
表面 先后 经 历 了 逆 压 梯度 和 顺 压 梯 度 的 作用 ， 流 体 
在 60% Cu 后 再 次 受到 逆 压 梯度 影响 ， 直 到 83% Ca。 
位 置 摩擦 力 小 于 零 ， 边 界 层 在 此 位 置 发 生 了 分 离 ， 
97% Ca 后 摩擦 力 变 为 正 值 且 急 剧 增加 ， 意 味 着 边界 
层 内 的 流动 状态 由 层 流 转换 为 湛 流 ， 并 再 附 到 壁面 
上 ,83%~97% Ci 之 间 是 回流 区 , 其 长 度 即 为 时 均 的 
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分 离 泡 长 度 。 另 外 ，0 ~ 5% Ci 之 间 摩 擦 力 为 负 值 ， 
说 明 由 于 较 大 的 正 攻 角 ， 使 得 前 缘 吸力 面 侧 出 现 了 
一 个 小 的 分 离 泡 。 其 他 四 个 工 况 中 摩擦 力 系数 的 变 
化 趋势 与 +7.8° 攻 角 类 似 ， 但 零 攻 角 和 人 负 攻 角 工 况 均 
没有 出 现 前 缘分 离 的 现象 。 随 着 攻 角 的 减 小 ， 分 离 
点 向 下 游 移 动 ， 而 再 附 点 位 置 基 本 没有 变化 ， 均 在 
97% Cu 附近， 因此 可 以 判断 ， 分 离 泡 尺度 也 随 着 攻 
角 的 减 小 而 缩短 。 

针对 -10>、0° 和 +7.8° 攻 角 ， 沿 叶 村 吸力 面 取 五 个 
位 置 ， 分 别 为 5% Cj,、30% Cj、70% Cj.、80% Ca 
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图 6 吸力 面 切 应 力 时 空 分 布 图 


Fig.6 S-T diagram of wall shear stress on the suction surface 
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图 7 吸力 面 时 均 摩 擦 力 系数 分 布 


Fig.7 The distribution of time-averaged friction coefficient 
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和 90% Ci， 将 时 间 平 均 和 展 癌 平均 后 的 二 维 速度 矢 
量 分 解 为 切 向 速度 和 法 向 速度 ， 法 线 上 的 切 向 速度 
在 图 8 中 给 出 。 可 以 看 到 ， 在 前 缘 5% Cu 处， 三 
种 攻 角 工 况 速度 剖面 相差 较 大 ， 且 攻 角 越 大 ， 法 线 
上 的 切 向 速度 越 大 ;， 从 壁面 附近 切 向 速度 的 局 部 放 
大 图 (图 9(a)) 中 可 以 看 出 ，+7.8° 攻 角 条 件 下 ， 靠 
近 壁 面 很 注 的 区 域内 出 现 了 回流 现象 ， 而 零 攻 角 和 
负 攻 角 工 况 均 没有 回流 发 生 。 随 着 法 线 位 置 向 下 游 
移动 ， 不 同 攻 角 下 速度 剖面 的 差异 逐渐 减 小 ， 直 到 
70% Ca 处， 各 工 况 下 的 ww 分 布 基本 一 致 。 由 图 7 可 


5%C 30%C。 TO%Cs 80%C 90%C 
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图 8 吸力 面 不 同位 置 的 切 向 速度 分 布 (一 代表 -10° 攻 角 ， 
代表 0° 攻 角 ，… 代 表 +7.8° 攻 角 ) 
Fig.8 Tangential velocity profiles at $ cross-sections on the 


suction Slide (一 : i= -10°, ---: 1= 0°, *…: = +7.8°) 
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0.006 


全 0.003 


(b) 90% Cu。 


图 9 5%Ca 和 90%Cax 处 近 壁 面 切 向 速度 (一 代表 -10° 攻 角 ， 
--- 代 表 0° 攻 角 ，… 代 表 +7.8° 攻 角 
Fig.9 Tangential velocity profiles at 5% Cand 90% Ci near the 


suction Slide (一 : i= -10°, ---: 1= 0°, *…: = +7.8°) 


— 


知 , 从 80% Cu 开始 各 攻 角 工 况 的 边界 层 将 依次 发 生 


分 离 ， 因 此 速度 剖面 又 表现 出 差异 。90% CGC, 处 边界 
层 内 速度 分 布 的 区 别 在 图 9(b) 中 更 清晰 地 展示 出 来 ， 
-10? 攻 角 条 件 下 ， 回 流 区 中 负 的 切 向 速度 值 较 小 ， 边 
界 层 的 厚度 也 最 薄 。 
4 结论 

本 文采 用 自 编 程 三 维 可 压缩 的 大 涡 模 拟 方 法 ， 
详细 研究 了 低压 涡轮 吸力 面 边界 层 的 演化 过 程 ， 分 
析 了 来 流 攻 角 对 边界 层 分 离 再 附 的 影响 。 主 要 得 出 


以 下 


界 层 
本 文 发 展 的 大 涡 模 拟 方 法 具有 较 高 的 精度 ; 


离 泡 ， 边 界 层 发 生 分 离 后 ， 


结论 : 
(1) +7.8" 攻 角 下 叶 栅 表 
分 离 及 再 附 位 置 均 与 实验 数据 吻合 较 好 ， 


压力 系数 、 吸 力 面 边 
说 明 


(2) 吸力 面 尾 缘 在 某 些 时 刻 会 同时 存在 多 个 分 


依次 展现 出 A 涡 、 发 卡 涡 


等 旋涡 


的 裔 


多 态 ， 发 卡 涡 的 拉 伸 作用 导致 了 层 流 边 界 层 
溃 ， 最 终 转 换 为 清流 ， 并 再 附 到 壁面 上 ; 
《3) 当 来 流 从 正 攻 角 向 负 攻 角 变 化 时 ， 吸 力 面 


边界 层 分 离 点 向 下 游 移动 ， 再 附 点 基本 保持 不 变 ， 


分 离 


微弱 


泡 长 度 减 小 ; 
(4) 较 大 的 正 攻 角 条 件 下 ， 吸 力 面前 缘 会 出 现 
的 分 离 现 象 ， 吸 力 面 后 部 分 离 泡 厚度 也 较 大 。 
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